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ABSTRACT 

In 1977, Flowers and Ruderman described a perturbation that destabilises a purely dipolar 
magnetic field in a fluid star. They considered the effect of cutting the star in half along 
a plane containing the symmetry axis and rotating each half by 90° in opposite directions, 
which would cause the energy of the magnetic field in the exterior of the star to be greatly 
reduced, just as it happens with a pair of aligned magnets. We formally solve for the energy of 
the external magnetic field and check that it decreases monotonically along the entire rotation. 
We also describe the instability using perturbation theory, and show that it happens due to the 
work done by the interaction of the magnetic field with surface currents. Finally, we consider 
the stabilising effect of adding a toroidal field by studying the potential energy perturbation 
when the rotation is not done along a sharp cut, but with a continuous displacement field 
that switches the direction of rotation across a region of small but finite width. Using these 
results, we estimate the relative strengths of the toroidal and poloidal fields needed to make 
the star stable to this displacement and show that the energy of the toroidal field required for 
stabilisation is much smaller than the energy of the poloidal field. We also show that, contrary 
to a common argument, the Flowers-Ruderman instability cannot be applied many times in a 
row to reduce the external magnetic energy indefinitely. 

Key words: magnetic fields - (magnetohydrodynamics) MHD - stars: magnetic field. 



1 INTRODUCTION 

Large-scale magnetic fields are known to be present in a wide va- 
riety of stellar objects, meaning that in these stars the dipole com- 
ponent (together perhaps with some other low-order multipoles) is 
not much weaker than the rms surfac e field. Th e initial discovery 
of such fields was in Ap stars jB abcockl 1 1 9471) . Since then, they 
have been observed or inferred to exist in white dwarfs, neutron 
stars, upper-main-sequence stars, and in the central stars of plane- 
tary nebulae. These fields appear to be long-lived, since they do not 
evolve on a timescale a ccessible to observations (see, for instance, 
iBraithwaite et all20ld) . 

A common feature of all these objects is that, over most of 
their interior, they are stably stratified. White dwarfs and neutron 
stars have no significant convective region^ while upper-main- 
sequence stars only have a small convective core. Dynamo effects 
are therefore expected to be irrelevant in keeping the strength of 
the magnetic field constant. The maximum surface magnetic fluxes 
observed in all these objects are similar, $ ma x = 7ri? 2 i? max ~ 
10 27 ' 5 G cm 2 , where _B max is the highest surface dipole strength 
detected in each class of objects. These two features are consid- 
ered compelling arguments in favour of flux freezing during stellar 



1 Recently formed neutron stars are only convective for some seconds, and 
white dwarfs have a thin convective region near their surface. 



evolution, although magnetic flux is certain to be lost in supernova 
explosions, so it is plausible that the magnetic field is regenera ted 
somehow in a young neutron star dThompson & Duncan 1 1993h. 
The ratio of fluid to magnetic pressure is dReisenegger l2009h 



8ttP 



8n d GM 2 
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3 x 10° 



M 



(1) 



which is a very large number even for the most strongly magnetised 
stars, unless the internal magnetic fields are substantially larger 
than the surface magnetic fields. Since this ratio is so high, we 
do not expect these fields to significantly modify the hydrostatic 
structure of the stars (i.e. their internal density and temperature 
profiles). However, they can play a major role in their evolution, 
e.g., through the transport of ang ular momentum (see, for instance, 
iHeeer. Wooslev & SpruitfaoOOh . 

Even though these long-lived fields have been known to ex- 
ist for more than half a century, it has not been possible to find 
an analytic model for a field that is in a stable equilibrium. How- 
ever, stable configurations that do not seem to evolve on timescales 
comparable to the Alfv en crossing time have been fo und to exist via 
numerical calculations dBraithwaite & Spruit 2004]), where an ini- 
tially random field usually evolves into an approximately axisym- 
metric configuration that is a combination of toroidal and poloidal 
components of similar energies (for axisymmetric fields, toroidal 
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Figure 1. Displacement produced by the Flowers-Ruderman instability. The 
thick arrows are field lines from the dipole, the dashed line shows the plane 
along which the star is cut, and both halves are rotated in opposite directions 
as shown by the thin arrows. 



and poloidal refer to the azimuthal and meridional components of 
the field, respectively). 

The stability of purely pol oidal or purely toroidal fields has 
also been studied in the past. iTavler I dl973l) , using the energy 
method, suggested that every purely toroidal field is unstable on 
an Alfven timescale, indepen dent of the streng th of the field . 
iMarkev & Tavler 1 1 197311 19741) and independently IWright I ( fl973h 
discovered that purely poloidal fields with closed lines contained 
inside the star are affected by an instability that is very similar to 
the kink instabilities in a Z-pinch. 

A simple argument given by iFlowers & Ruderman I i 19771) 
shows that any purely poloidal field with field lines extending out- 
side the star should be unstable. If the initial configuration is such 
that the external field resembles a dipole, cutting the star in half and 
rotating each half by 90° in opposite directions (as shown in Fig. 
QJ would greatly reduce the dipole component of the field, lead- 
ing to a magnetic field with less energy. However, neither Flowers 
and Ruderman nor anyone else have given a formal proof of this 
argument. 

In the numerical simulations of iBraithwaite I d2009l) . instabil- 
ities related to the poloidal and toroidal components of the field 
are studied. Using the stable configurations found after simulating 
the evolution of random fields, Braithwaite used different ratios of 
poloidal to total energy of the magnetic field, Ep / E, and found the 
field to be stable for 0.056 < E P /E < 0.8. The field became un- 
stable for Ep/E > 0.8, with an m = 2 mode that seems to consist 
mostly of displacements in latitude of the fluid. For Ep /E > 0.9, 
modes with higher m became unstable, as would be expected since 
these modes have to overcome a higher resistance from the toroidal 
field. These modes resemble kink instabilities, as was mentioned 
before. 

The structure of this paper is the following: In S|2] we formally 
prove the Flowers-Ruderman instability for a pure dipole field. To 
do so, we explicitly calculate the energy of the external magnetic 
field as a function of the angle of rotation of each half of the star, 
and see that it is a monotonically decreasing function. In $3] we 
treat the problem using perturbation theory for a particular family 
of fields. In fj4] we estimate the stabilising effect of a toroidal field 
when the perturbation is not done with a sharp cut through the star, 
but rather with a displacement field that switches continuously from 
one direction of rotation to the other, over a thin but finite region. 
In ij5] we show that, contrary to what has been claimed previously 
in the literature, higher order multipoles cannot be achieved with 



successive cuts in different directions, and in f|6] we present the 
conclusions of our work. 



2 PROOF OF THE FLOWERS-RUDERMAN 

INSTABILITY BY AN EXACT EVALUATION OF THE 
ENERGY 

We consider the star to be non-rotating and completely surrounded 
by vacuum. Inside the star, the magnetic field is axisymmetric and 
purely poloidal, and outside the star the field is a pure dipole. If 
we completely ignore the effects of the magnetic field on the hy- 
drostatic structure of the star (an assumption we justify in ! j2.5l l. 
then the star should be perfectly spherical, and when the Flowers- 
Ruderman instability takes place, each half of the star rotates as a 
rigid solid. Since in stellar interiors the magnetic Reynolds number 
is much larger than 1, field lines will be dragged by the fluid with- 
out modifying the magnitude of the magnetic field at each point, 
and thus the internal magnetic energy of the star will not be modi- 
fied in the procesfl Therefore, we are only interested in the energy 
of the external magnetic field, and we now proceed to prove that 
this energy is in fact monotonically reduced by performing the dis- 
placement suggested by Flowers and Ruderman. 



2.1 Exterior energy of an arbitrary magnetic field 

To start, we must obtain the magnetic field outside the star, given 
the field on its surface. Because outside the star there are no cur- 
rents, we have V x B = and therefore B = V\£ r . Since 
V • B = 0, * must satisfy Laplace's equation, 



V 2 * = 0. 



(2) 



The solution to this equation in spherical coordinates that is of 
physical significance to us is 



(3) 



l — 1 m— — I 



where the term with / = is excluded since it corresponds to a 
magnetic monopole. The a; m are obtained by requiring the field to 
have a continuous normal component at the surface, 



dim = 



Y l * m (8,&(B r ) r=R dQ, 



(4) 



where B r = B ■ r is the radial component of the field. Now, 
the magnetic energy inside the star should not change, since the 
field only rotates while keeping its magnitude. However, the exte- 
rior field changes significantly. Thus, the variation of the magnetic 
energy can be obtained just by computing the variation outside of 
the star. The exterior magnetic energy is obtained from 



E = 



r 87T 



dV = 



(Vfr) 2 
8tt 



dV = 



V ■ (jgW) 
8tt 



dV (5) 



where V covers all space outside the star, and we used the fact 
that satisfies Laplace's equation. Using the divergence theorem, 



2 It is important to note that the plane that cuts the star cannot cross any 
field lines, otherwise, to rotate each half these field lines should be cut, and 
that is not possible. 
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the energy can be expressed as a surface integral, with a normal 
pointing into the star|j 

£= ILv[ ( v$) r=fi -dS. (6) 
8 71 " J s 

Since we consider the star to be perfectly spherical, 
V* • dS = -R 2 (B r )r=RdQ, and consequently 

" 8tt 



E = 



4n S7r i4vr 



Q (7) 



where in the last step, we used the fact that <]/ is real and set it equal 
to its conjugate. Replacing the expression for ^, as given by eq. 
l|3)> we get the following result: 



I m 
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~8tv 



Y — 
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ai m {l + 1) 
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(8) 



where we used eq. Q to express the integral in terms of the aim. If 
we define c !m as 



aimjl + 1) 
R>+ 2 



YC m {6,cj>){B r 



sdtl, 



the energy of the external magnetic field is 

p i? 3 |c !m | 2 



8tt ^ f + 1 

lm 



(9) 



(10) 



Another useful result is the radial component of the field at the sur- 
face expressed in terms of the c; m coefficients. This can be obtained 
from the radial component of V^, or by inverting eq. {5), and it is 
equal to 



(B r )r=R = 2^ c tmYtr. 



(ID 



2.2 Proof that the final energy is less than the initial one 

The results contained in eq. j 10b can be used to prove that the en- 
ergy is effectively reduced when one half of the star is rotated with 
respect to the other. To do so, let us define a quantity T as 

f {B r f r=R dQ. 



T 



8tt 



(12) 



This quantity will be conserved when the star is cut in half and 
rotated. So, using the superscripts i and / to denote initial and final 
states, T 1 = T f . If we use the result of eq. il lb to express the terms 
of (B r )l =R = (B r B*) r=R , we get 



T 



El 

8tt 
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By rewriting T 1 = T f , we obtain 

E l C i"i| 2 = E l C '"i| 2 - 



(13) 



(14) 



(15) 



3 The term corresponding to the surface at infinity vanishes. This is because 
the term in that decreases most slowly with r goes like r _2 ,so(>I'V'I 1 )- 
d S goes like r -3 and vanishes in the limit r — > oo. 



If the initial external field is a dipole field, the only nonzero cj m is 
c\ . Considering this, 



E i c '" 



l c iol 



From here, using eq. d 1 Ob we get that 



E t < 



8tt 
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R 3 \c[ \ 

8tt 2 



Ei. 



(16) 



(17) 



Thus, the final state will have less or equal energy than the initial 
one. The equality would hold if and only if the c\ m are equal to zero 
when I ^ 1, which is not the case since the severe discontinuity 
that is produced cannot be resolved into an expansion of spherical 
harmonics with a finite number of terms. 

It is important to note, however, that we only proved that the 
magnetic energy of any final state after cutting the star and rotating 
it is less than the initial energy of the dipole field. We have yet 
to prove that the energy is monotonically decreasing for the entire 
rotation. So, up to this point, we could expect the minimum energy 
to be present at some intermediate point in the rotation, and not 
after the rotation has been completed. 

2.3 Proof that the energy decreases monotonically 

The initial field outside of the star is that of a dipole, so 

(B r ) r =R = B cose, (18) 

where Bo is the strength of the field at the poles. This can be written 
in the form of eq. Q lb as 



/ 4.7T 

(B r ) r =R — Boy _ 3"^' 10 ' 



(19) 



so from eq. l[9j, the only nonzero c; m is cio = Bo\/47r/3, and 
using eq. dlOb , the external magnetic energy of this field can be 
evaluated as 



E c 



BpR 3 
12 ' 



(20) 



Now, we rotate each half of the star by an angle 6 in opposite 
directions, 



(B r 



B cos[0'(0,<t>,G)] x>0 

B cos[o'(e,4>, -e)] x<q 



(21) 



where 0' corresponds to the polar angle in a spherical coordinate 
system that rotates together with each half of the star, as illustrated 
in Fig. [2] It is fairly straightforward to show that 



cos[0 (0, (j>, O)] — cos 6 cos O — sin 6 sin sin 0, 
so that the radial component of the field at the surface is 

J Bo (cos 9 cos — sin 8 sin </>sin 0) x > 
I Bo (cos 6 cos + sin 6 sin <fi sin 0) x < 

For this field, the c; m as given by eq. l|9} are 



(22) 



(23) 



cim = cos 0Bo / cos 9Y lm d Q 

J4tt 

r 37r/2 



+ B sin 



sin 9 sin (f>Y lm d (f> d I 



Jtt/2 



(24) 
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Figure 2. Relationship between the coordinates 9 and 0' used in i]2.3l Ini- 
tially we have a spherical coordinate system in which the polar angle 9 of a 
point in space corresponds to the angle formed with respect to the z axis. A 
new axis z' is defined by performing a right handed rotation by an angle © 
with respect to the z axis, and the polar angle 9' of a point is defined as the 
angle formed with respect to the z' axis. 



and due to the symmetries of Yi m , it can be shown that 



2V7r/3B cose i = l,m = 
cim = ^ 2wi m Bo sinO I and m even 
otherwise 



(25) 



where 



3tt/2 



wi m = I dO d^>sin 9 siti(f>Yi m . (26) 

/o Jtt/2 

Using these values for the c; m combined with eq. d 10b . the external 
magnetic energy is obtained as a function of B, 



E = cos 2 QE C + sin 9 



R 3 B 2 
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Defining 
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TV 

lm 



Wlr, 



l + l 



the energy can be rewritten in a compact form as 

E = E c [1 + sin 2 Q(A - 1)1 . 



(28) 



(29) 



Since the complete rotation is obtained with B = tv/2, the energy 
is a monotonic function of B. If A > 1, the energy will be an 
increasing function, but if A < 1, it will be a decreasing function. 
However, from the results of $\2.2\ we already know that the final 
energy is smaller than the initial one, thus A < 1 and the energy 
decreases monotonically along the entire rotation. Therefore, the 
Flowers-Ruderman instability is present in the case of the purely 
dipolar field. 

Even though we already proved the existence of the instability, 
an estimate of A is called for. To obtain this estimate, we consider 



the quantity 



* = ; E E 



TV ^ — I' + l 

l'—2 m=—l 
l' even rn even 



(30) 



which tends asymptotically to A as I increases. Also, due to the 
conservation of T described in E12.2I and the initial and final values 
of the c; m , which can be obtained by setting B equal to zero or 90° 
respectively in eq. {25}, it can be seen that the wi m must satisfy 



\V)lm\ = IT- 



(3D 



Using this, we can obtain lower and upper bounds on the value of 

A, 



Ai < A < Ai + 



6 



tt(7 + 3) 



;-e e 



\ 



l'=0 rn=—l' 
l' even rn even 



2 



(32) 



and since both limits tend asymptotically to A as I increases, this 
in principle can be used to evaluate A to an arbitrary precision. By 
using / = 100, one finds that 



0.5463 < A < 0.5466. 



(33) 



The final energy of the system is Et = AE C , so t he process nearly 
halves the energy of the external magnetic field. iRoberts I J 198 lh 
gave a value of A — 0.577, which is slightly higher than ours, but 
he did not describe exactly how he computed that number. 



2.4 Including higher order multipoles in the initial 
configuration 

Extending the previous analysis to consider the superposition of 
multipole components with the dipole component turns out to be 
a much more complicated problem. However, it is possible to get 
a simple answer if we only consider axisymmetric fields, and if 
we only care about the energy difference between the final and the 
initial state. 

Consider an axisymmetric field that consists of a dipole and 
some higher order multipole, of fixed order I > 1, 



* = 



(f + l)- 



■l+l 



-Y- 

1 10' 



(34) 



B =V#, 



where the superscript i is used to denote that these coefficients rep 
resent the initial state. After the rotation, the field might not be ax 
isymmetric, and it will be defined by some set of coefficients cf m 
From eq. d 1 Ofc . the final energy of the system is 

p3 |„f |2 

Et 



8i7 ^ 



\ c lm\ 
l+l 



(35) 



It can be shown that the final configuration has no dipole com- 
ponent, by proving that each multipole by itself does not produce 
a dipole component after the rotation. To do this, we consider one 
half of the star to be rotated by 180° instead of both halves being 
turned by 90° in opposite directions as was done in the previous 
section. The choice of how we perform the rotation is irrelevant 
to our calculation, since the energy associated with each multipole 
cannot depend on the particular spherical coordinate system chosen 
to perform the spherical harmonic expansion. However, the true 
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physical process requires both halves to rotate in opposite direc- 
tions due to the conservation of angular momentum. 

In particular, it is evident that the final configuration has no 
dipole component when I is even, since in these cases Y^ is sym- 
metric with respect to the equator, and a rotation by 180° of one 
half of the star will leave the field just as it was at the beginning. 
For the case of odd I, the dipole components c[ m that would be 
produced by that particular multipole after the rotation can be ex- 
pressed from eq. ^ as 

c im = A tm / P lm (cos<9) • P fo (cos6Osin0d0, (36) 
Jo 

where the Pi m are associated Legendre polynomials 
dArfken & Weber I l2005h and Ai m contains the result of inte- 
grating over <j>, the normalising factors of the spherical harmonics, 
and an additional factor that relates to the strength of the multipole. 
Due to the orthogonality condition for the associated Legendre 
polynomials, c[ is equal to zero, and since Pi±i(cos6) ~ sin# 
(except for a numerical constant), we have 



f 

Cl±l 



P fo (cos6>)sin 8dd 



=A t±1 J ^ P io (x)y/l-x*dx = 0, 



(37) 



where the equality to zero is due to the fact that the P;o ( x ) are °dd 
functions of x for odd I. Thus, the dipole component is zero after 
the rotation, so eq. {35} can be written as 
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Due to the conservation of T shown in eq. d!5t . 
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(40) 



where Eo and Em correspond to the initial energies of the dipole 
and the multipole components, 



E D = 



R 



3 'clol 2 



Em = 



8tt 2 ' 8vr / + i ' 

and the initial energy of the system is simply Ei = Eo 



(41) 



Em- 



Combining eqs. d38t and d40t a sufficient condition for the final 
energy to be smaller than the initial one (i.e. Ei < Ei) can be 
obtained, 



-Eo 



l + l 



Em < Eo + E^ 



Eo 

Em 



>l 



(42) 



(43) 



This, however, does not ensure that the initial state is unstable, since 
the energy might not necessarily be monotonically reduced along 
the rotation, but it shows that a state with lower energy can exist, 
even when a significant fraction of the energy is contained in a 
higher order multipole. This is particularly true for the quadrupole, 
in which case eq. J4-3 b shows that independent of the strength of the 
quadrupole with respect to the dipole, the final state still has lower 
energy. 

The condition given in eq. d43t might look very restrictive for 
multipoles with a very high /, but this is only because it is a suffi- 
cient but not a necessary condition. 



2.5 Validity of the assumption of sphericity 

Although the overall change in the structure of the star due to the 
presence of the magnetic field is small, it is not obvious that the en- 
ergy associated with this perturbation is smaller than that released 
by the instability just discussed, and can therefore be ignored. To 
check if this could in a way stabilise the star against the Flowers- 
Ruderman instability, we consider the same displacement, where 
each half is rotated with respect to the other. As the magnetic en- 
ergy inside the star and the internal energy of the fluid will remain 
the same, we need only to consider the changes in the external mag- 
netic energy, and the gravitational potential energy. 

The change in the structure of the star due to the presence of 
the magnetic field is related to the plasma /3 parameter (which in 
our case is /3 2> 1). For instance, the displacement of the surface of 
the star compared to the perfectly spherical one is of order R. 
The change in the external field between both situations is due to 
the small perturbation of the surface, which scales like f) , so, if 
we use B to denote the field in the case where the perturbation 
to sphericity is ignored, and B + AB to denote the case where 
the perturbation to sphericity is considered, we should have that 
AB/B ~ The energy associated with the external field, if 
we were to consider the change in the structure of the star, would 
only be a correction of order /3^ 1 E C , so, it can be safely ignored. 

We require a slightly more intricate argument to show that the 
change in the gravitational potential energy during the rotation can 
be ignored. Consider the star as composed of four different mass 
distributions, by considering each half separately, and furthermore, 
by taking each half as composed of a spherical distribution of mass 
plus a small perturbation due to the presence of the magnetic field. 
The spherical distribution would have a density p, while the pertur- 
bation Ap would be of order /3 _1 p. The total gravitational poten- 
tial energy of the system can then be decomposed in several terms, 
four corresponding to the individual gravitational energy of each 
mass distribution, and twelve more related to the interactions of 
each mass distribution with the other three. Of course, when per- 
forming the rotation of one half of the star with respect to the other, 
the gravitational energy of each of the four elements remains the 
same, and the only possible change is due to the interactions be- 
tween them, which can be treated case by case, as follows: 

(i) First, we take the potential energy due to the interaction 
among the spherical mass distributions. When one half rotates with 
respect to the other, the system is equivalent to the initial one, so 
this does not produce a change in the potential energy of the sys- 
tem. 

(ii) Next, we consider the interaction between a spherical mass 
distribution and the perturbation of a spherical mass distribution. 
In this case, a rotation of one half leaves the system composed of 
these two parts in a state that is equivalent to the initial one, and dif- 
fers only by a rotation of the whole system. Thus, the gravitational 
energy in this case remains the same. 

(iii) Finally, we take the interaction between the two perturba- 
tions to the spherical mass distributions. In this case, after the ro- 
tation, the system is generally not equivalent to the initial one. The 
mass associated with each distribution is ~ /3 _1 M, so the potential 
energy associated with their interaction is of order 



,GM 2 ! B 2 GM 2 

= /3 



R 



R 



and, since P ~ GM 2 /R 4 , 
E g ~ f3~ 1 B 2 R 3 ~ /3 _1 £ c 



(44) 



(45) 
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so the change in the gravitational potential can in fact be ig- 
nored compared with the change in the external magnetic energy. 
Of course, an exact solid body rotation of each half is not re- 
ally expected to happen when the star is not spherical, since that 
would produce some "rough edges". But, since the perturbation of 
sphericity can be ignored in that case, we expect a more reasonable 
displacement that closely resembles the rotation of each half with 
respect to the other to exist and reduce the energy even further. 



3 PROOF OF THE INSTABILITY USING 
PERTURBATION THEORY 

Using MHD perturbation theory, we should also be able to prove 
the existence of the Flowers-Ruderman instability. This proof how- 
ever could not be as complete as the one given in ij2] since perturba- 
tion theory can only be used to see if the system is unstable against 
small displacements, and thus, we cannot prove with this approach 
that the energy decreases monotonically along the entire rotation. 
Nevertheless, we now provide a proof of the instability for a par- 
ticular family of fields using perturbation theory, since the results 
obtained in doing so will be useful in the next Section, where we 
will require the result for the energy perturbation in terms of the 
angle by which both halves of the star are rotated. 

In £13-11 we prove that the volume contribution to the potential 
energy perturbation is equal to zero. In S13.2I we calculate the con- 
tribution to the potential energy perturbation due to surface currents 
and show that the final result directly relates with the energy given 
in eq. {29}. 



3.1 Contribution to the potential energy perturbation inside 
the star 



Using the energy principle of iBernstein et al. 1 i 1958h . the stability 
of a system perturbed by a displacement field $ is given by the sign 
of the potential energy perturbation, which can be written as a sum 
of hydrostatic and magnetic terms 

SW =5Whyd + Wmag, 



SWhyd = 



[riP(v ■ £) 2 + (£ ■ vp)(v-0 

-(£ ■ V$)(V • P £) + p£ ■ V5<f>] d V 



£ ■ (Sj x B + j x SB) d V 



(46) 



where V now denotes the volume of the star, P is the fluid pressure, 
p is the mass density, $ is the gravitational potential, Ti is defined 
as 



Fi 



d\nP 



(47) 



and j is the current density, which in a static configuration can be 
written as 

4.7T 

— j = V x B. (48) 

c 

The magnetic field and current perturbations are given by 

47T 

5B = V x (£ x B), — Sj = VxSB. (49) 

c 




Figure 3. Cylindrical coordinate system used in ij3.ll The coordinates of 
this system are ro, •& and z', and the cylinder is oriented in such a way that 
z' = x and that for 1? = one has ■oj = z. Also, the basis vectors of 
this coordinate system satisfy the relation z! = va X ij, which defines the 
direction of increasing 



If SW < 0, then the resulting configuration will be unstable. We 
ignore the effects of the magnetic field on the structure of the star, 
so P, p and <E> are spherically symmetric. 

The displacement field for the case of the Flowers-Ruderman 
instability is taken to be 



Qrx x r = —Or (cos 8 cos <f>cj> + sin cj>9) x > 
-Qrx x r — Or(cos 6 cos (f><f> + sin (f>6) x < 



(50) 



with |0 1 <C 1. This displacement field has no radial component 
and is incompressible, so <5Whyd = 0. 

It can also be proved that for this displacement field, 
fSWmag = in the bulk of the star, independent of the form of 
the magnetic field. This is to be expected, since the magnetic field 
will be simply displaced, and its energy should not change. To do 
this, we consider a cylindrical coordinate system with coordinates 
vj for the cylindrical radial coordinate, # for the azimuthal angle 
and z' , as illustrated in Fig. [3] In each half of the star, the displace- 
ment field can be written as 

£ = ±©c?i9, (51) 

and in this system, a general magnetic field can be written as 

B = B„m + B^ + B z ,z', (52) 

where B^, B$ and B z i are all 27r-periodic functions of $ and must 
be such that V • B = 0, i.e. 



1 d 



(ujB^) + 



1 dB# dB z 



0. 



w dvj vo d-d dz' 

Using this, the integrand in eq. J46t can be written as 

S. i6j xB+jx5B) = ±-—, 

where F is a 27r-periodic function of # given by 



F = B Z 



f dB z 



dBj 
1 dz' 



+ B Z 



fdB^ 



8B 6 



(53) 



(54) 



(55) 



Since the range of integration for $ in eq. ( 146b is from to 27r, and 
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the integration is done over the derivative of a 27r-periodic function, 
<5Wmag = in the bulk of the star. 

However, we already saw that there is an effective variation 
of the energy when performing this perturbation, and thus, we are 
not taking into account all the work that is done on the fluid. This 
large scale displacement produces surface currents in two different 
regions, and these should be responsible for the work done: 

(i) Along the surface of the sphere. Since the exterior field sat- 
isfies Laplace's equation, and its boundary conditions only require 
the normal component of B to be continuous, it is unlikely that a 
large-scale displacement that affects the surface of the star will not 
produce a discontinuity of the tangential component of B in some 
areas. Also, even before the rotation is made, surface currents will 
be present if the field is discontinous along the surface of the star. 
Thus, surface currents are an important element for perturbations 
that affect the surface. 

(ii) Along the plane that cuts the star. The discontinuity pro- 
duced by the rotation will produce a current sheet along this plane. 

From these two effects, only the first is really relevant to the energy 
of the star. The second effect is not, because £, B, and j are parallel 
to that surface, and thus no work is done on the fluid. 



to SW due to these terms can be written afl 

' 2 



SW S c = 



[£ • (SK x B + K xSB) 



(61) 



However, due to the discontinuity of the tangential components of 
B and SB across the boundary, the choice for these two vectors 
is somewhat ambiguous, a consequence of the unphysical nature 
of surface currents. This can be avoided by considering only per- 
turbations that are parallel to the surface, so £ r (R, 8, <fr) = 0, in 
which case only the radial components of B and SB contribute to 
the previous expression, which reduces to 



5W SC = - 

87T 



B r £-(SB cxt -SB int )dQ 



SB r £-(B cxt -B int )dn 



(62) 



Here, it is not necessary to distinguish between the interior and 
exterior values of B r and SB r because these must be continuous. 
The primary difficulty in this expression is the term <5B oxt . How- 
ever, this integral can be explicitly computed in terms of the 5ci m . 
To do this, we write £ = ^cft + ^gO, B — B r r + BgO, and using 
integration by parts it can be seen that 

?3 |r„ |2 



8tt 



E 



\Scin 



l + l 



(63) 



3.2 Contribution to the potential energy perturbation due to 
surface currents 

If a discontinuity of the 6 and cj> components of the magnetic field 
exists at the surface, a s urface current K will be produced (see, for 
instance. I Jackson 1998), 



47T 

— K = f X (B c 

c 



B in 



(56) 



where Bi nt and B cxt are the fields immediately inside and imme- 
diately outside the surface of the star respectively. If the field is 
perturbed by a displacement £, then Bj n t changes to first order in 
£ by 5Bi n t — V x (£ x 2?i n t). This change will modify the bound- 
ary conditions on the stellar surface, giving rise to a perturbation of 
the exterior magnetic field 



<5B cxt = V<5#, 
with 



(57) 



v- R 5c lm 



(1 + 1) 



Ylr, 



Sci m = / j£,(M)(tfBr)r=*dn, 



(58) 



where 8B r — SB ■ r is the perturbation of the radial magnetic field 
component. Almost certainly, this will give rise to a perturbation of 
the surface current: 



47T 

— SK = r x (<5Bext - SB lnt ). 



(59) 



Now, by replacing j by j + 5(r — R)K in eq. ( 146b and performing 
the radial integral for the term with the unperturbed surface current 
and the one with the surface current perturbation, the contribution 



so this term is always positive, and thus does not drive the instabil- 
ity. 

Now we consider the perturbation field given by eq. d50fc and 
a dipolar magnetic field as given by eq. d 1 8 b - In this case, the Sci m 
arc 



2B &w u 




I, m even 
otherwise 



(64) 



where Bo is the strength of the field at the poles. Using this, to- 
gether with eq. ( 163b . we get 



f B r £ ■ SB cxt dn = AE e e 2 , 



(65) 



where A is given by eq. ( 1281 ) and E a is given by eq. l !20t . This gives 
us one of the terms of SW SC (as shown in eq. J6>2b> ■ The other terms 
can be evaluated directly by using the displacement field given by 
eq. d50b and a magnetic field given by eq. i l 1 St . yielding the result 



8tt 



+ / SB r £. ■ (B c 



Bi„ t )dfi 



(66) 



Considering eqs. {62}, J65l l and \66\ . the potential energy perturba- 
tion due to the surface currents is found to be 



5W BC = (A- l)E e e, 

which agrees with eq. J29t up to order O 2 , as expected. 



(67) 



4 Considering this, SWmug now consists of a volume integral and a surface 
integral: 



SW, 



mag---, £-(Sj xB + jx6B)dV 

* Jv 



2 J4w 



f [£ • (SK x B + K x SB)] r=R dfl 

JAtt 



(60) 
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4 STABILISING EFFECT OF A TOROIDAL FIELD 

When a toroidal field is added, magnetic field lines will pass 
through the plane that cuts the star in half in the Flowers-Ruderman 
instability. Thus, if a sharp cut is done, magnetic field lines would 
be cut, which is not possible. Because of this, an arbitrarily weak 
toroidal field is enough to stabilise the star against the sharp cut, 
but if the cut is done smoothly, toroidal field lines will not be cut, 
but instead will be severely twisted. To study this, we consider the 
effects of performing the cut of the star smoothly across a region 
of finite width 2eR (as shown in Fig.[4j. As e increases, this bend- 
ing will be less pronounced, and thus the stabilising effect of the 
toroidal field will be reduced. Under some reasonable assumptions, 
we use perturbation theory to obtain a ratio between the energy 
of the poloidal field and the total energy of the magnetic field for 
which the field becomes stable to this displa cement. Thi s value can 
be compared with the values obtained by Brai thwaite 1 J2009I) for 
which the field becomes unstable in MHD simulations. 

To do this, we consider a displacement field of the form 



-Qorx x f x < —eR 
£ = { Q(x)rx x r \x\< eR 
Qorx x r x > eR 



(68) 



where Q(x) is a continuous, odd function in the interval \x\ < eR 
that satisfies 



Q(±eR) = ± Q , 



dQ 

dx o 



: 0. 



(69) 



The condition imposed on the derivative is to avoid discontinuities 
in SB along the boundary, which would in turn produce surface 
currents. Similar to the displacement field used before for the sharp 
cut, £ has no r component and satisfies V ■ £ = 0, so there will be 
no hydrostatic contribution to SW as can be seen from eq. J46b . 

The potential energy perturbation for this displacement field 
can be split into several terms, including a term that involves sur- 
face currents, 8W SC - This contribution to the potential energy per- 
turbation involves surface integrals of an infinite number of spheri- 
cal harmonics, and the fact that the displacement field is defined in 
terms of Cartesian coordinates adds great complexity in trying to 
evaluate SW SC - Because of this, we consider that the smooth tran- 
sition is done in a thin region relative to the radius of the star, so 
e< 1, and we assume that 5W S c does not change significantly with 
respect to the value obtained for the sharp cu{f|. In any case, we ex- 
pect SW SC to increase as e increases, since in this case the dipole 
component of the external magnetic field will not be reduced as 
much as in the case of the sharp cut. 



4.1 Structure of the toroidal field 

In an axisymmetric hydromagnetic equilibrium, the magnetic field 
and its associated currents cannot produce forces in the azimuthal 
direction, because the pressure and gravity forces would not be able 
to balance them. This implies that the toroidal magnetic field is 



5 We do not expect the external magnetic field to be significantly different 
on the surface of the star for the region \x\ > eR, so the contribution to 
SWsc on this region should not change significantly. Also, the area of the 
surface in the region \ x\ < eR is small compared to the rest of the surface in 
which the integral for 8W B c is done, so even if there are significant changes 
there, we do not expect them to significantly modify the work done on the 
whole surface. 




Figure 4. Smooth rotation of each half of the star due to the inclusion of a 
toroidal field. The thick arrows are field lines from the dipole, and the thin 
lines indicate the rotation of each half of the star. The dashed lines enclose 
the region of width 2eR where the displacement field switches the direction 
of rotation continuously, as is shown by the thin arrows there. 



contained in a particular region of the star that is restricted by the 
topology of the poloidal field, as will be shown in this section. 

The most general axially symmetric magnetic field can 
be decomposed into a toroidal and a poloidal part, each of 
which is determined by a scalar function (see, for instance, 
IChandr asekhar & Prendergast |[l95r3) . 

B = Bt + Bp, B T =/3V0, Bp = Va x V</> (70) 

whert0/3 = P{r, 9) and a = a(r, 6). Since Bp _L Va and B T -L 
V/3, a and /3 are constant along poloidal and toroidal field lines 
respectively. The force per unit volume exerted by the magnetic 
field is given by 



cF M =j x B 



=j P x Bp + j P x B T + j' t x Bp +j T x B T , 



(71) 



where j p and j T are the currents related to the poloidal and 
toroidal fields, 



4.7T 



-JP 



■In 



V x Bp, — j T = V x B T - (72) 

c ' c 

Since j P and Bt are both toroidal fields, their cross product is 
zero. Also, j p x Bp and j T x Bt are both the cross products 
of a toroidal and a poloidal field, so they have no <j> component. 
The remaining term, j T x Bp is the cross product of two poloidal 
fields, so it points in the <j> direction, and the <j> component of the 
magnetic force is given by 



cF M <t,4> = j T x Bp. 



(73) 



If the configuration is in axisymmetric equilibrium, Fm^ must van- 
ish since there is no possible way for the fluid to counteract this 
magnetic force, so j T || Bp. The current j T can be calculated as 



4%j T = V x B T = V/3 x V</>, 



(74) 



so the condition that j T |] Bp is equivalent to the condition (V/3 X 
V</>) || (Va x W(f>), and since V/3 and Va have no <fi component, 
V/3 || Va which means that a and f3 can be written as functions 
of each ot her, a = a(/3) or = /3(a). This res ult was originally 
derived by Chandrasek har & Prendergast I dl956h . 

6 In this subsection we use /3 to represent the function that defines the 
toroidal field, instead of its usual meaning as a ratio of fluid to magnetic 
pressure. 
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Now, consider a poloidal field line that closes outside a star 
that is in equilibrium. Along this field line, a is constant, and thus, 
ft — /3(a) is also a constant. However, outside the star the field can- 
not have a toroidal component, so /3 = in that region, and since 
P is constant along the field line, the toroidal field vanishes every- 
where along a poloidal field line that closes outside the star. Be- 
cause of this, the toroidal field must be contained in regions where 
the poloidal field lines are closed within the star. This is not a new 
result, since this condition for axisymmetri c hydromag netic equi- 
libria has been known for many years (e.g. Roberts Ill98ll) , 

In Fig.|7]we plot the field lines of the particular poloidal com- 
ponent of the magnetic field we use in this section (this field is 
described in greater detail in £|4.3.2b . For this poloidal field, the 
toroidal field is enclosed in a torus-like region, so over most of the 
star the field is purely poloidal. 







f J 













4.2 Cylinder approximation and toroidal fields 

As a simple approximation to the region of transition (|x| < eR), 
we will consider it as a cylinder of height 2eR and radius R, and 
use the same cylindrical coordinate system used in §3.l\ 

The displacement field of eq. J68t in this region can be written 

as 



£ = 9(z')ot?. 



(75) 



We consider the perturbation in the potential energy of a toroidal 
field due to this displacement. Since the height of the cylinder is 
small relative to the radius of the star, we approximate the toroidal 
field as 



S T = &0,tf)z' 



(76) 



where b(zu, ■£)) is a 27r-periodic function that is odd in $ (i.e. 
b(zu, —•&) — —b(zzs, 1?)). Using this, the potential energy pertur- 
bation in this region due solely to the toroidal field can be obtained 
from eq. ( 1461 ) . 



SWt = - — 

07T 



dz' 



n2 d ( db\ 



dzu 



ZU B 



' dz' 2 



(77) 



The first term vanishes after integration over #, and since we de- 
mand that the derivative of <d vanishes on the boundary, the second 
term can be rewritten after integration by parts as 



SWt 



1 

8^ 



dz' 



~dz' 



J 3,2 

dzu zu b . 



(78) 



From this, it can be seen immediately that SWt > 0, so, as ex- 
pected, the toroidal field opposes this displacement. 

We now need to specify a model for both Q(z') and b(w, 
We choose our function O(z') as 



O(z') =Oo sin 



SWt 



HZ 

2^R 
32eR 



(79) 



'o Jo 

This function Q(z') is odd and satisfies the required conditions 
mentioned in eq. J69l >. With this particular displacement field, 
SWt oc e _1 , so, as mentioned before, if the region where the dis- 
placement field switches direction is very thin, the magnetic energy 
will increase significantly, and thus an infinitely weak toroidal field 
is enough to stabilise the star against a sharp cut. 



Figure 5. Model used for the toroidal field. The vertical line is the symme- 
try axis, and the field is contained in a circular torus of radius fiR that is 
tangent to the equator of the star, as shown by the shaded region in the fig- 
ure. Also depicted in the figure are the coordinates p and ip used to describe 
the magnitude of the field. In the cylinder approximation the torus is treated 
as two independent cylindrical regions. 



Since the toroidal field is confined within the poloidal field 
lines that are closed inside the star (as is shown in §4.11 , we con- 
sider the toroidal field to be contained in a torus of internal radius 
pR. This approximation is adequate for the particular poloidal field 
we will use, as can be seen in Fig. [7] In the cylinder approximation, 
we consider this torus as two cylindrical regions of radius pR that 
are centred at (zu, #) = (R(l — /J.), ±7r/2), as illustrated in Fig. [5] 
In each of these regions, the strength of the field will depend on the 
distance to the centre, so we switch to coordinates (p, ip) centred 
on one of these circles in which we have b = b(p) (as is shown in 
Fig. |5} . The corresponding SWt can be solved in these coordinates 
as 



SWt = 



IQeR 



dp dp pb (p)d (p,<p), 



(80) 



where d 2 (p, ip) 



[R(l - p)f - 2pR(l - p) cosip is the 



distance to the origin. As a model for b(p), we use 



b(p) = r/Bo cos' 



(81) 



where Bo is the maximum strength of the poloidal field on the sur- 
face, and r/Bo is the maximum strength of the toroidal field. The 
square on the cosine assures SB to be continuous along the surface 
where the toroidal field vanishes. Using this model for the field, 
SWt results in 



SWt 

Hp) 



3E c Ql rj 2 h(n) 



64tt 2 
-■R 2 (6n 2 



32)p(I - 
+ (9tt 4 



2p) 

- 77tt 2 + 192)//, 



(82) 



where E c is the initial energy of the exterior magnetic field as given 
by eq. j20b . 

It can be seen from eq. d80t that the detailed geometry of the 
toroidal field is not so relevant, specially if the toroidal field is con- 
tained in a region far away from the centre of the star. In the latter 
case, d(p,ip) ~ R(l — p), and the integral will involve only the 
square of the magnitude of the magnetic field times an area ele- 
ment. Because of this, SWt should be closely related to the energy 
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Figure 6. Plot of the function g(p) from eq. i84t 



of the toroidal field, rather than depend on its detailed geometry. 
The energy of the toroidal field can be computed as 

£ t = ^(37r 2 -16)?7V(l-M), (83) 
which can be used to rewrite eq. ( 182b as 



(3tt 2 - 16) /i 2 (l -/x)' 

A plot of the function <?(a0 is shown in Fig. [6] where it can be seen 
that at a fixed energy of the toroidal field, the stabilizing effect will 
be stronger if the field is contained farther away from the center of 
the star. Also, since g(fi) only varies by a factor ~ 2, it is the energy 
of the toroidal field what sets the order of magnitude of SWt, while 
its geometry plays a secondary role. 

4.3 Effect of poloidal fields for the smooth rotation 

4.3.1 Cross term in SW 

When a poloidal field is added, a cross term appears in SW that 
involves both the poloidal and toroidal components of the magnetic 
field. This term has the form 



SWcioss = 



\ j dV£- [Sj T x B P + j T x SB T 



(85) 



+Sj P x B T + j P x 8B T ] ■ 



Considering only the parity of the functions involved, it can be 
shown that the integrand in (JWa-oss is an odd function of x, and 
since the integral is over the interval —eR<x< eR, integration 
over x will immediately give zero as a final result, so 

SW CIOSS = 0. (86) 



4.3.2 Purely poloidal contribution to SW 

For the poloidal field, we will consider a configuration of the form 



B = Vq x V4>, a = f(r) sin 



(87) 



On the surface of the star, the radial component for these fields 
is 2f(R)R~~ 2 cos(#), and thus, outside the star all these fields 
are pure dipoles. This model for the internal field covers a wide 
range of axisymmetric configurations; this includes the constant 
field B — Bqz in which f(r) = £>or 2 /2 and the fields used by 




Figure 7. Field lines for the particular poloidal field configuration we use 
in this section. The dashed line represents the surface of the star, the solid 
lines are the field lines, and the region where poloidal field lines are closed 
inside the star is marked in grey. This field is of the form given by eq. {87), 
with f(r) given by eq. (88). The region where poloidal field lines are closed 
inside the star has a shape similar to a torus, and the toroidal field must be 
contained there. Figure adapted with permission from Akgiin et al. (paper 
in preparation). 



iBraithwaite I d2007h to study the stability of purely poloidal fields 
in rotating stars. 

Our choice for the function f(r) is equivalent to that of Akgiin 
et al. (paper in preparation), but normalised so at the poles the 
strength is Bo- This poloidal field is of the form of eq. {87}, with 



f(r) = 



35B 
16 



_ 6/ 3 
" 5R 2 + 7R 1 



(88) 



This field is completely continuous across the surface of the star, 
so there are no surface currents present in the equilibrium config- 
uration. Also, it satisfies \j \ = at the surface, which is expected 
from the fact that the matter density goes to zero there. A plot of 
this field is shown in Fig. [7] 

Using the cylinder approximation, it is difficult to treat the 
contribution to SW due only to the poloidal field. It is also diffi- 
cult to treat the problem in spherical coordinates, since the regions 
of integration involved are non-trivial. However, for our particular 
choice of the poloidal field, the purely poloidal contribution to SW 
can be obtained exactly using Cartesian coordinates. The displace- 
ment field used here is of the form of eq. {68}, and the function 
Q(x) is the same as that of eq. {79} with z' replaced by x. 

With all this, the potential energy perturbation due solely to 
the poloidal field can be computed from eq. {46}, which was done 
using the software Maximal]. The result is a finite polynomial in e, 
to lowest order 



SW P 



(23tt 2 



8192 
-0.15E e Qle. 



330) B 2 R 3 e 2 e. 



(69tt 2 



2048 



"Ve 2 , 



(89) 



This contribution is negative, but it is not as important as that of 



7 http://maxima.sourceforge.net 
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Figure 8. Integrand for SWp, — £ ■ (Sj P X Bp + j P X 5Bp)/2, in the 
plane given by x = e/2 using e = 1/10, with the line where the integrand 
is equal to zero plotted on top. The regions that contribute negatively to 
SWp (black in the figure) correspond closely to the region where poloidal 
field lines are closed inside the star, meanwhile the regions that contribute 
positively to SWp (white in the figure) lie near to the symmetry axis. 



SW SC (from eq. d67t it can be seen that 5W BC ~ -0.45£ c e 2 ). Ini- 
tially we expected the poloidal field to perform a stabilising effect, 
since this displacement would tend to twist field lines that are near 
to the symmetry axis. However, the region where the poloidal field 
lines are closed within the star turns out to be highly unstable to 
this displacement, as can be seen in Fig. [8] It can be seen that the 
contribution to the potential energy perturbation is positive along 
the axis of symmetry, and the region where it is negative encloses 
the field lines that are closed inside the star. We believe the positive 
contribution to be caused by the twisting of field lines, and the neg- 
ative contribution to be due to the displacement of closed poloidal 
field lines, which resembles a kink mode as described by Markey & 
Tayler (1973) and Wright (1973). In a kink mode, closed poloidal 
field loops inside the star are displaced with respect to each other 
mostly in a direction perpendicular to the neutral line, but still, the 
displacement parallel to the neutral line is fundamental. Perhaps 
adding such component to our analysis could cause an important 
decrease in SWp, but this would significantly complicate our anal- 
ysis, and it is unlikely that we can find an analytic solution for SWp 
under those circumstances. 



4.4 Total potential energy perturbation 

To obtain the total energy perturbation, we add all the contributions 
obtained so far, given by eqs. (|67j>, (84} and {89]>, 



5W = SWsc + SWt + SWp 



■E e @ 2 



-(1-A) + 



E T 990 - 69tt 2 



E c 8vre 



2048 



(90) 
(91) 



If SW — 0, then the system is marginally stable, and for that case, 
solving for Et /E c in terms of n and e results in 



Et 

~eZ 



87T6 

9(t*) 



(990 - 69tt 2 
2048 



(92) 



which is an increasing function of both fi and e in the range that 
these variables cover (0 < (X < 0.5, < e < 1). Choosing n 
and e, we obtain a lower bound on the strength of the toroidal field 



needed to stabilise the star against a smooth rotation done over a 
region of width 2eR. However, /i is not completely arbitrary, since 
in equilibrium, the toroidal field must be contained by the field lines 
that are closed inside the star. A reasonable value for fj, (for the 
poloidal field chosen) is fj, — 0.2, and slight variations from this 
value do not modify the result significantly. 

The choice for e must be a small value in order for our ap- 
proximations to be valid. We choose in a rather arbitrary way the 
value e = 1/3 for which the region in which the displacement field 
switches the direction of rotation has a width of 2R/3. Under these 
conditions the assumption that 8W SC is similar to the one calcu- 
lated for the straight cut in i j3.2l is very doubtful, so this choice of e 
certainly serves as an upper bound. 

Using n — 0.2, e = l/3 and A — 0.546 (which is taken from 
eq. J33t ) to evaluate Et/E c in eq. {92}, one obtains Et/E c ~ 
0.27. 



4.5 Comparing the poloidal and toroidal energy of the 
magnetic held 

In order to compare this result with that of Braithwaite (2009), we 
must see what it means in terms of the energies of the toroidal and 
poloidal fields (for the latter, including the volume outside the star). 
The energy of the poloidal field can be evaluated as 
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With this, the ratio of poloidal to total energy is 

Ep Ep 1 

~~E ~ E T + E P ~ 11 E T ' 
70 E c 

For the value Et/E b ~ 0.27 obtained in the previous section, 
this ratio is very close to unity, Ep /E ~ 0.96. This tells us that a 
toroidal field with an energy much smaller than the poloidal field 
is enough to stabilise the star against this perturbation. This can be 
co mpared with the inst ability that could be seen in the simulations 
bv lBraithwaitel |2009) for a ratio of E P /E = 0.8. As this per- 
turbation happens with a much stronger toroidal field, it seems to 
indicate that the perturbation we are studying is not the dominant 
one, since other instabilities are present for the poloidal field even 
when the toroidal field is strong enough to stabilise it against the 
one we have studied. 



5 EFFECT OF SUCCESSIVE CUTS 



Some authors, for instance, iRoberts I dl98lh and 
iBraithwaite & Nordlund I |2006), have stated that a second 
cut in a direction perpendicular to the first one would produce 
a configuration which resembles an octupole. Additionally, they 
expected that this process could be repeated ad infinitum, with cuts 
in different directions, to produce a configuration that resembles 
an arbitrary multipole. Since the external energy of the field is 
expected to be reduced as higher order multipoles are achieved, in 
principle they expected the process to proceed naturally, reaching 
very high order multi poles, for which the external energy of 
the field is negligible. IRoberts I d 198 lh computed the energies 
associated with the external magnetic field for several of these 
multipoles, showing that the energy of each multipole was in fact 
smaller than that of the previous one. Even though our result differs 
slightly from the value obtained by Roberts for the quadrupole, we 
expect this tendency to be true. 
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Figure 9. Performing two cuts to the dipole. The star is seen from the top 
of the symmetry axis, and the plus and minus signs indicate magnetic lines 
coming out and into the star respectively. At the top the wrong picture is 
shown, where an octupole is produced, while the bottom shows that actually 
the second cut leaves the star in the dipole configuration. 




Figure 10. Performing several cuts to the dipole in order to obtain an oc- 
tupole configuration. 



However, this process breaks down already at the second cut, 
as is shown in Fig. [9] where it is seen that the second cut leaves 
the star again in the quadrupole configuration, instead of producing 
an octupol^fl It is possible to obtain the octupole with a sequence 
of several cuts, as shown in Fig. [10] but through this process the 
energy does not seem to be monotonically reduced, since interme- 
diate steps have important contributions from a dipole component, 
so the star cannot actually follow this sequence of displacements. 
Even though we only showed a particular sequence of cuts that 
transforms the quadrupole into an octupole, clearly this cannot be 
achieved by a single cut, and it does not seem plausible that this can 
be done without increasing the external magnetic energy at some 
point, so we expect this mechanism to affect only the initial ax- 
isymmetric field. 



8 Note that the terms quadrupole and octupole used here do not refer to a 
pure quadrupole or octupole, but to a configuration in which the primary 
component is that of a quadrupole or an octupole respectively. In effect, 
these configurations have contributions from an infinite number of multi- 
poles. 



6 CONCLUSIONS 

iFlowers & Ruderman I dl977h presented an argument that shows 
how purely poloidal fields in stars are unstable. If the external field 
is similar to a dipole, one could cut the star in half and rotate each 
piece in opposite directions, leading to a configuration in which the 
external field resembles a quadrupole, and thus, the energy of the 
external magnetic field should be significantly reduced. Although 
the Flowers-Ruderman instability is widely accepted, no formal 
proof had been given that shows both that the external magnetic 
energy is reduced when the rotation of each half is completed and 
that the energy decreases monotonically along the entire process. 

In this work, we presented a formal proof of this mechanism 
for the case in which the initial field outside the star is that of a 
point dipole, by computing the energy of the external field along 
the entire rotation. We showed that the external magnetic energy 
decreases monotonically, having a final value of approximately 
0.55-Ec, where E c is the initial external magnetic energy as given 
by eq. i20l . 

We also studied the Flowers-Ruderman instability using per- 
turbation theory, in which case we had to consider the effects of 
surface currents in order for the instability to appear. These effects 
are not unique to the Flowers-Ruderman instability, and should be 
considered for any displacement that modifies the magnetic field 
on the surface. The result obtained for the potential energy pertur- 
bation of the star was found to be consistent with the exact value of 
the energy previously found. 

We then studied how a toroidal field could stabilise the star 
against the Flowers-Ruderman instability. Since a sharp cut through 
the star would split toroidal field lines, the rotation has to be car- 
ried out with a continuous displacement field that switches the ori- 
entation of rotation across a thin region. For a specific model, it 
was found that the configuration was stable against the Flowers- 
Ruderman instability for a ratio of poloidal magnetic energy to to- 
tal magnetic energy of Ep/E < 0.96. Using MHD simulations, 
Braithwaite] d2009l) had shown that when the ratio Ep /E was be- 
low 0.8, the instabilities driven by the poloidal field were sup- 
pressed, but if the ratio was just above 0.8, the field was found 
to be unstable with an m = 2 mode that does not resemble the 
Flowers-Ruderman instability. Because of this, we conclude that 
the Flowers-Ruderman instability is not the dominant one. 

However, the critical value for the ratio Ep /E we obtained is 
not a very accurate measure of the stability of the system against 
the perturbation we treated, since the contribution to the potential 
energy perturbation comes from different sources which do not nec- 
essarily scale with the energy of the poloidal field. If a configura- 
tion only has a small region with closed poloidal field lines, then 
the kinks that produce the negative contribution to the energy de- 
scribed in SI4.3.2I will not be important. In a similar way, the mag- 
netic field can be very strong inside the star, but weak at the sur- 
face, so the energy released by the Flowers-Ruderman instability 
can change by orders of magnitude between configurations with 
similar poloidal magnetic energy. In the case when the effect due to 
the closed poloidal fields is not important, a good indicator should 
be the ratio of external magnetic energy to toroidal magnetic en- 
ergy, which for our case results in Et/E ~ 0.27 as the critical 
value, which is computed from eq. i92\ with the values previously 
obtained for /i and r\. 
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